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Abstract 

In this series of lectures the basic ideas of the 1/rriQ expansion in QCD (rag is the mass 
of a heavy quark) are outlined. Applications to exclusive and inclusive decays are given. 
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1 Introduction 



With the precise formulation of the 1/ttiq expansion in QCD using effective field theory and 
operator product expansion heavy quark physics has been based on model independent ground, 
which allows to reduce the uncertainties due to the QCD bound state problem drastically. From 
the point of view of weak interactions the main interest in processes with heavy quarks lies in 
the exploration and the test of the CKM sector of the standard model (SM) describing the 
masses and the mixing of quarks. From the experimental side the running experiments as well 
as the ones planned in the near future constitute a large effort to explore this part of the SM, 
which is not yet tested with an accuracy comparable with the one of the coupling of the Z to 
the fermions. 

The main theoretical progress in the description of systems involving a single heavy quark 
is based on the infinite mass limit of QCD 0, in which two additional symmetries appear that 
are not present in full QCD ||. This limit may be regarded as the leading term of a 1/tuq 
expansion and a systematic approximation to full QCD may be constructed using the methods 
of effective field theory, the so-called Heavy Quark Effective Theory (HQET) Q. The heavy 
mass limit and applications of HQET have been extensively studied, and the development of 
the field is documented in more or less extensive reviews [[|. 

The corrections to the heavy mass limit are characterized by two small parameters, namely 
the strong coupling constant, taken at the scale of the heavy quark « s (mg) and the ratio A/rriQ 
of the scale of the light degrees of freedom A and the heavy quark mass. While the first kind of 
corrections may be calculated perturbatively in terms of Feynman diagrams, the second type 
needs additional non-perturbative input encoded in the matrix elements of higher dimensional 
operators. 

However, the non-perturbative input is to be taken at a small scale, where the symme- 
tries of HQET may be applied and hence the non-perturbative physics is constrained by these 
symmetries. In particular, for heavy to heavy decays these additional symmetries restrict the 
number of independent form factors to only a single one, of which the absolute normalization 
at a specific kinematic point is fixed by the symmetries. 

By combining the method of the 1/wq expansion with the short distance expansion one 
may obtain a heavy mass expansion also for inclusive decay rates ||-[[7|,||-|l3]. The heavy 
quark mass sets a scale that is large compared to Aqcd, and one may use a similar setup as 
in deep inelastic scattering for the description of inclusive decays. In this way one may not 
only study total rates, but also differential distributions such as the lepton energy spectra in 
inclusive semileptonic decays. 

This series of three lectures tries to summarize the basic ideas of the heavy mass expansion. 
In the first lecture the formulation of HQET as an effective field theory is described and the 
additional symmetries of heavy quarks are introduced. The second and third lecture deal with 
applications to exclusive and inclusive decays. 



2 Effective Field Theory 



Effective field theories jnj have become a widely used tool in modern elementary particle 
physics. An effective theory treatment is convenient if the problem under consideration involves 
very disparate mass scales such that the physics that is to be described happens at much lower 
energies than the scale set by some heavy particles in the theory. In such a case it is useful to 
switch to an effective theory in which the heavy degrees of freedom do not appear explicitly; 
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they only reappear in the effective theory as higher dimensional operators, which are multiplied 
by coupling constants with negative mass dimension. The scale of the coupling constants is set 
by the large mass and thus these contributions are small, if the scale of the physics described 
with the help of the effective theory is small compared to this large mass. 

An effective theory is always valid only in a limited region of scales, a natural cut-off is 
given by the mass of the particle which has been removed by switching from the full to the 
effective theory. As mentioned above an effective theory involves interactions which would 
lead to a non-renormalizable theory, if one would consider the theory to all orders in these 
higher dimensional operators. However, working to a definite order of the expansion in inverse 
powers of the large scale one does not face any problem concerning renormalization. Starting 
from the renormalizable dimension-4 piece of the effective theory Lagrangian we may use its 
renormalization group properties to study the cut-off dependence of the effective theory, which 
is determined by the short distance properties of the effective theory. 

Applying effective theory methods corresponds to an expansion of the Greens functions of 
the full theory in inverse powers of the large mass scale; such an expansion is only possible up to 
logarithmic dependences on this large scale. These logarithms may be accessed using a properly 
constructed effective theory, where these logarithms correspond to the renormalization group 
logarithms of the cut off. In this way one may even achieve a resummation of the logarithmic 
terms using renormalization group methods in the effective theory. 

In the case at hand the large scale is the mass tuq of the heavy quark, and the leading term of 
this expansion is the static limit. This effective theory (Heavy Quark Effective Theory, HQET) 
is a a powerful tool, allowing for numerous purely QCD based calculations. Renormalization in 
this effective theory implies a factorization theorem for the Greens function of full QCD, which 
means that to any order in the I/tuq expansion one may factorize the short distance physics 
from the long distance effects. Explicitly this means for a Greens function Gj u u calculated in 
full QCD that one may rewrite it as 

G fu u{p 1 ---p n ,m Q ,fi = m Q ) = Y,[—) Z (j \m Q ,fi)G ( ^ f (p 1 ---p n ,fi). (1) 

j \ m QJ 

where /xq is the renormalization point of the full QCD function. Here the constants Z^\rriQ, //) 
depend on the factorization scale \x and on logarithms of itiq\ it contains all the short distance 
effects, which may be calculated perturbatively. Using the renormalization group of the effective 
theory, one may perform a systematic resummation of logarithmic dependences on the heavy 
quark mass. 

The Greens functions G^jj are calculated in HQET and do not depend on the heavy mass 
any more; they contain the long distance effects which are not calcualable via perturbation 
theory. However, as we shall see below, they are constrained by heavy quark symmetries. 
Consequently one gains complete control over the mass dependence by switching from full 
QCD to HQET, the effective theory obtained in the heavy mass limit. 

In this section we shall discuss the formulation of this effective theory. First we shall review 
and compare some of the different possibilities to formulate the infinite mass limit. In this limit 
new symmetries appear which are the key to various model independent statements concerning 
weak decay matrix elements; these symmetries are reviewed in paragraph 3.1. Finally, we shall 
consider the systematic approach to the calculation of corrections to the infinite mass limit. 
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2.1 Infinite Quark Mass Limit as an Effective Field Theory 



This issue has been discussed repeatedly and various formulations of the infinite mass limit 
are available. Of course, as far as physical quantities are concerned, all approaches yield the 
same result. However, for some special applications one approach may be more convenient than 
another. 

The equivalence of all different approaches is ensured by a theorem well established in the 
field of effective theories. It has been shown that in an effective theory involving a field <j) one 
may perform redefinitions of the fields , such that 

0^P(0,«90) (2) 

where P is an arbitrary polynomial function |]I5| , |T6| . 

Such a redefinition will not change the S'-matrix, although the Lagrangian (and Greens 
function) expressed in the redefined fields may look completely different. In this sense the 
different formulations of HQET are equivalent. 

In the following we shall consider two formulations of HQET. The first one nicely exhibits 
the fact that HQET is an effective theory in the sense that one integrates out a heavy degree of 
freedom and performs an expansion of the remaining action functional in the large mass scale, 
which for the case at hand is the heavy quark mass. The process of integrating out the heavy 
degree of freedom may be performed explicitly as a Gaussian functional integral, and one may 
construct a formulation of HQET be expanding the result in powers of I/itlq. In this way one 
obtains a 1/ itlq expansion of both the heavy quark field as well as for the Lagrangian. 

The second formulation is based on the standard way of separating "upper" and "lower" 
components of the spinor fields by performing a sequence of Foldy Wouthuysen transformations. 
These transformations lead to an expansion in 1/rriQ for the heavy quark field and for the 
Lagrangian, which is different for each of these quantities from the other formulation. However, 
if one calculates a physical quantity, both approaches will yield the same answer, because the 
1/rriQ expansion has to be unique for an observable. 



Integrating out heavy degrees of freedom 

One may obtain a formulation of the heavy mass limit by integrating out heavy degrees of 
freedom from the functional integral of QCD Greens functions [17]. This integration may in 
fact be done explicitly, since for the case at hand it amounts to a Gaussian functional integration. 
Formulating the heavy mass limit in this way clearly exhibits that it corresponds to an effective 
theory in the usual sense. We start from the generating functional of the QCD Greens functions 

Z(r], fj, A) = J [dQ][dQ)[d(j)x] exp jzS + iS x + i J d 4 x (fjQ + Qrj + (f>\\) j , (3) 

where = q, A ^ denotes the light degrees of freedom (light quarks q and gluons A^) with the 
action S\, while S denotes the piece of the action for the heavy quark Q including its coupling 
to the gluons 

S = J d A xQ(ilp-m Q )Q. (4) 

where 

is the covariant derivative of QCD. We have introduced source terms r\ for the heavy quark and 
A for the light degrees of freedom. 
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We shall consider hadrons containing a single heavy quark, and we assume that this heavy 
hadron moves with a velocity v 

v = Madron ? y 2 = ^ Vq > o (6) 
m hadron 

This velocity vector may be used to split the heavy quark field Q into an "upper" component 
and a "lower" one x 

4> v = ^(l + ^)Q, t<j) v = (t), (7) 

Xv = \(l-i>)Q, tXv = -X, (8) 

and to define a decomposition of the covariant derivative into a "longitudinal" and a "trans- 
verse" (_L) part 

D fl = v ll (vD) + D^, D^ = {g^- Vlx v v )D\ , j) = 0. (9) 

Using (0-0) the action (^) of the heavy quark field takes the form 



S = (Fx 



(j){i(v-D)-m Q }(p-x{i(v-D)+mQ}x + (j>ip X + xip <P 



(10) 



The heavy quark in the meson is very close to being on shell, and thus the space time dependence 
of the heavy quark field is mainly that of a free particle moving with velocity v. This suggests 
a reparametrization of the fields by removing the space time dependence of a solution of the 
free Dirac equation. We shall choose the "particle-type" parametrization corresponding to the 
"positive energy solution" of the Dirac equation 

0, = e- im ^h v , X v = e~ m ^H v , (11) 

such that the space time dependence of the remaining fields h v and H v is determined by the 
residual momentum k = p — mqv, which is due to binding effects of the heavy quark inside the 
heavy hadron, and which is a "small" quantity of order Aq C d- 
Expressed in these fields the action of the heavy quark becomes 

S = J d 4 x [h v i{v ■ D)h v - H v {i(v ■ D) + 2m Q }H v + h v iIp L H v + H v Hp L h v ] , (12) 

The term containing the sources is also rewritten in terms of the fields h v and H v 

j d A x (fjip + tjjr]) = J d A x {p v h v + h v p v + R V H V + H V R V ), (13) 

where p v and R v are now source terms for the upper component field h v and the lower com- 
ponenent part H v respectively. 

In terms of the new variables the generating functional reads 

Z(p v ,p v ,R v ,R v ,X) = J [dh v ][dh v ][dH v ][dH v ][d<i>x] (14) 
exp liS + S\ + i J d 4 x (p v h v + h v p v + R V H V + H V R V + 0aA)| , 



where the action S for the heavy quark is given in eq . (|12|) . 

From ( |12D it is obvious that the heavy degree of freedom is the lower component field H v , 
since it has a mass term 2mg, while the upper component field h v is a massless field describing 
the static heavy quark. In the heavy mass limit only the Greens functions involving the field 
h v have to be calculated, and hence we integrate over H v in the functional integral fll4]) with 
the sources of the lower component field R v and R v set to zero. This can be done explicitly, 
since it is a Gaussian integration 



Z(p v , p v , A) = / [dh v ] [dh v ] [d\)A exp liS + S\ + i / d 4 x (p^h^ + h^p^ + A A) 



(15) 



where now the action functional for the heavy quark becomes a non-local object 

1 



S = id x 



hti(vD)K-htp- U— y — — 



This Gaussian integration corresponds to the replacement 

/ 1 



(16) 



2rriQ + iv D 



ip±h v 



(17) 



for the lower compponent field. Furthermore, the Gaussian integration yields a determinant 
A. In the full theory one may also perform this Gaussian integration, and the determinant 
obtained contains all the closed loops of heavy quarks. After renormalization of the full theory 
their contribution starts at order 1/m 2 with an Uehling potential like term. In the effective 
theory one may take the determinant A to be a constant, if the terms of order 1 /trig and higher 
comming from the closed heavy quark loops are included by matching to the full theory. Since 
we shall discuss only the leading term of the I/tuq expansion in this section, we may drop the 
determinant in what follows. 

The non-locality of the action functional is connected to the large scale set by the heavy 
quark mass, and the non-local terms may be expanded in terms of an infinite series of local 
operators, which come with increasing powers of I/itlq. In the context of a field theory this 
corresponds to a short distance expansion and hence these operators have to be renormalized. 
The tree level relations may be read off from the geometric series expansion of the non-local 
term in (jlB]). In this way we obtain the expansion of the field and the Lagrangian 



Q(x) 



c 




(18) 



h v (ivD)h v + h v ip 

1 

2m 



(19) 



h v {ivD)h v + -^h v {ip ± ) 2 ih v + (— ) h v {ip ± ){-ivD)(ip ± )h v + 



A Greens function with an operator insertion is treated in a similar way; the heavy quark 
fields entering the inserted operator are dealt with in the same way. The net effect of this is 
that the heavy quark fields in the operator insertion are replaced by the expansion 
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Foldy Wouthuysen Transformation 

A second way of formulating the heavy mass limit proceed along the well known steps performed 
in deriving the non-relativistic limit of the Dirac equation ||18|| . The reasoning used here is 
motivated by quantum mechanics; as a first step one rewrites the equation of motion in a 
Hamiltonian form 

ivdQ = HQ, H = f{p L + m Q + gvA) (20) 

where A is the gluon field. Note that in the rest frame we have v = (1, 0, 0, 0), and thus (|20|) 
takes the usual form, since vd = do and v A = Aq. 

In general the Hamiltonian couples the upper and the lower component of the heavy quark 
field Q, projected out by P± = (1 ± f). The Foldy Wouthuysen transformation is a transfor- 
mation of the form 



Q^Q' = exp(iF)Q H -> H' = exp(zF) [H - ivd) exp(-iF) 



(21) 



where F is a hermitian matrix in the space of the Dirac spinors, such that exp(iF) is unitary 
in that space. This requirement is motivated by quantum mechanics where the spinor is in- 
terpreted as a wave function, but this interpretation becomes meaningless once we switch to a 
field theory. 

The transformation ( 21 ) is designed such that the resulting Hamiltonian H' does not couple 
upper and lower components of the field Q any more. The generator F of this transformation 
may be expanded in powers of l/mg, from which one may construct a l/mg expansion of the 
Hamiltonian and the transformed fields. Removing the mass term of the Hamiltonian by a 
phase redefinition as in (O), one obtains for the fields and the Lagrangian 



Q(x) 



£ = K 



-imqv-x 



iv ■ D — 



hJx) 



(22) 



2m, 



-Pi + 



Ami 



-fllv -D + p ± v Dp ± - l -v ■ Dpi 



(23) 



We note that the leading order term as well as the terms of order 1 / mQ are identical in the two 
approaches. Differences start to appear at order l/mg, which are terms involving a factor which 
would vanish by leading order equations of motion. In the Fouldy Wouthuysen formulation the 
Lagrangian does not contain such terms, while these terms appear in the approach of integrating 
out the lower components of the heavy quark field. However, subleading terms of a physical 
matrix element will consist of local contributions originating from the expansion of the field Q 
as well as of non-local pieces involving time-ordered products of the leading order currents with 
the subleading terms of the Lagrangian (see below). If the time-ordered products are taken 
with a term that would vanish by a naive application of the leading-order equation of motion, 
these will lead to a contact terms, i.e. effectively to a local contribution. In this way the terms 
in the first approach rearrange in such a way that the final result for an observable quantity 
is the same in both cases. For some practical applications the Foldy Wouthuysen approach 
has an advantage that all the time-ordered product with the Lagrangian are truly non-local 
contributions, in other words, none of the contributions will lead to a contact term. 
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2.2 Corrections to the Heavy Mass Limit 
Tree Level Considerations 

Corrections to the infinite mass limit may be considered in a systematic way. Starting from the 
tree level expressions given in the last sections one may use the expansion of the Lagrangian and 
the fields as given in section 2.1. to construct the 1/toq expansion of full QCD matrix elements. 
In doing this it will not matter which representation of HQET (e.g. the one that is obtained 
from integrating out the heavy quark or the one constructed from the Foldy Wouthuysen 
transformation), since the matrix elements have to have a unique I/tuq expansion. 

As an example we shall consider a matrix element of a current qTQ mediating a transition 
between a heavy meson and some arbitrary state \A). The full QCD Lagrangian £ and the 
fields Q are expanded in terms of a power series in 1/mg in the way described in section 2.1, 
and the matrix element under consideration up to order 1/mg takes the form: 

(A\qTQ\M(v)) = (A\qTh v \H(v)) (24) 
— (A\qTP_iph v \H(v))-i [ d i x(A\T{L 1 (x)qTh v }\H(v)) + C(l/m 2 ) 

n Q J 



2m 



where L\ are the first-order corrections to the Lagrangian as given in ([II]) or (p23[) . Furthermore, 
\M(v)) is the state of the heavy meson in full QCD, including all its mass dependence, while 
\H(v)) is the corresponding state in the infinite mass limit. 

Expression fl2~4] ) displays the generic structure of the higher-order corrections as they appear 
in any HQET calculation. There will be local contributions coming from the expansion of the 
full QCD field; these may be interpreted as the corrections to the currents. The non-local 
contributions, i.e. the time-ordered products, are the corresponding corrections to the states 
and thus in the r.h.s. of (0) only the states of the infinite-mass limit appear. If one switches 
to another representation of HQET, one reshuffles terms from the fields into the Lagrangian; in 
this way the Lagrangian picks up operators which are proportional to the equations of motion. 

As an example for the kinds of matrix elements appearing in subleading orders of the 1 / tuq 
expansion we consider the mass of a heavy hadron. In the infinite mass limit this mass is given 
in terms of the quark mass plus some "binding energy" A. The corrections are of order I/tuq 
and are given by the matrix elements of the leading correction term of the Lagrangian. One 
obtains |19 

m H = m Q 1 1 + — + — ^ (Ai + d H X 2 ) + 0{l/m%) ) (25) 

\ m Q Zm Q J 

where dn — 3 for the 0~ and dn = — 1 for the 1~ meson. Up to this order no non-local 
terms appear; such terms show up the first time at order l/m?Q. The parameters A, Ai and A2 
correspond to matrix elements involving higher order terms that appear in the effective theory 
Lagrangian 

(0|gjjjg isK\H{y)) 

(Q\q l5 h v \H{v)) 1 ) 

(H(v)\h v (tD) 2 h v \H(v)) 
A i - KT7 



2M 



. (H{v)\h v a^iDHD»h v \H(v)} 

As " 2M H (28) 
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where the normalization of the states is chosen to be {H (v)\h v h v \H (v)) = 2Mh = 2(mg + A). 
These parameters may be interpreted as the binding energy of the heavy meson in the infinite 
mass limit (A), the expectation value of the kinetic energy of the heavy quark (Ai) and its 
energy due to the chromomagnetic moment of the heavy quark (A2) inside the heavy meson. 
The latter two parameters play an important role since they parametrize the non-perturbative 
input needed in the subleading order of the 1/mg expansion. 

The only parameter which is easy to access is A2, since it is related to the mass splitting 
between H(v ) and H*(v, e). From the S-meson system we obtain 

A 2 (m b ) = ±(M H . - M H ) = 0.12 GeV 2 ; (29) 

from the charm system the same value is obtained. This shows that indeed the spin-symmetry 
partners are degenerate in the infinite mass limit and the splitting between them scales as 

The other parameters appearing in (p5|) are not simply related to the hadron spectrum. 
Using the pole mass for txiq in (|^), QCD sum rules yield for a value of A = 570 ± 70 MeV 
More problematic is the parameter Ai; from its definition one is led to assume Ai < 0; a more 
restrictive inequality 

- Ai > 3A 2 (30) 
has been derived in a quantum mechanical framework in E| and using heavy-flavour sum rules 



21 1 . Furthermore, there exists also a sum rule estimate [22] for this parameter: 



Ai = -0.52 ±0.12 GeV 2 . (31) 

This value is compatible with the bounds; however, it is unexpectedly large since it corresponds 
to a rms-momentum of the heavy quark inside the meson of 



(p 2 ) ~ 720 MeV (32) 

which is large compared to the naive guess of — Ai ~ (Aqcd) 2 This is the reason why also 
smaller values of Ai have been used in the literature. 

Recently there has been an attempt [^] to extract A and Ai from the shape of the lepton 
energy spectrum in inclusive semileptonic B decays (see sections 4.3 and 4.4). The values 
obtained from this analysis are A = 0.39 ± 0.11 GeV and — Ai = 0.19 ± 0.10 GeV 2 , where the 
MS definition of the mass has been used. The uncertainties quoted are only the la statistical 
ones; the systematical uncertainties of this approach are difficult to estimate. 

Beyond Tree Level 

Going beyond tree level will induce corrections of order a"(mg), n = 1,.... These may be 
calculated in terms of Feynman diagrams which may be evaluated using the Feynman rules of 
HQET. Only two Feynmal rules are modified compared to full QCD: 



full QCD HQET 

% i 

p — rriQ + ie vk + ie 

Heavy quark gluon vertex ig^^T 11 — > igv^T 01 



% % 

Propagator of the heavy quark — > — , p = mv + k 

p — rriQ + ie v k + ie 
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For the sake of clarity we shall stick to our example of a heavy light current considered above. 
To leading order in the 1/wiq expansion one may evaluate the radiative corrections to such a 
matrix element using the above Feynman rules and finds a divergent result with a divergence 
related to the short distance behaviour. Since HQET is an effective theory, the machinery of 
effective theory guarantees the factorization of long distance effects from the short distance 
ones, which are related to the large mass mg. Neglecting I/tuq corrections, this factorization 
takes the form 

(A\qTQ\M(v)} = Z (^j (A\qTh v \H(v))l + 0(l/m Q ) (33) 

From Feynman rule calculation one obtains the perturbative expansion of the renormalization 
constant Z which generically looks like 

Z {^j = a 00 (34) 
+ a n ^ s ln^^j\ +a 10 a s 

+ a 22 (a s In (^-^j J + a 2i«s (a> s In j + a 2o«^ 

+ a 33 (a s In ( ^ jj + a 32 a s ^a s \n ^^jj + a 31 a 2 s ( a s In jj 
+ a 30 a 3 s H 

where ot s = g 2 /(An). 

This factorization theorem corresponds to the statement that the ultraviolet divergencies 
in the effective theory have to match the logarithmic mass dependences of full QCD. The 
factorization scale /i is an arbitrary parameter, and the physical quantity (A\qTQ\M(v)) does 
not depend on this parameter. However, calculating the matrix element of this operator in the 
effective theory and studying its ultraviolet behaviour allows us to access the mass dependence 
of the matrix element (A\qTQ\M(v )). 

The ultraviolet behaviour of the effective theory is investigated by the renormalization group 
equations. Differentiating ( |3~3"D with respect to the factorization scale fi yields the renormaliza- 
tion group equation 

iL, {z(j) wWm)l} = (35) 

from which we may obtain an equation which determines the change of the coefficient Z when 
the scale is changed 

/ d +7 M)zM=0 (36) 



\d\iafi 

u^) = -^—H(A\qrh v \H(v)%). 

dmjj, 

The quantity 7j is called the anomalous dimension of the operator J = qTh v which is universal 
for all matrix elements of this operator, since it is connected with the short distance behaviour 
of the insertion of J. 

Eq.(|36D describes the renormalization group scaling in the effective theory. It allows to shift 
logarithms of the large mass scale from the matrix element of J into the coefficient Z: If the 
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matrix element is renormalized at the large scale ttiq the logarithms of the type lnmg will 
apear in the matrix element of J while the coefficient Z at this scale will simply be 

Z(l) = a 00 + a w a s (m Q ) + a 20 ^(m Q ) + a 30 a£(rag) H (37) 

The renormalization group equation (^) allows to lower the renormalization point from mq to 
fi; the matrix element renormalized at ji will not contain any logarithms of mg any more, they 
will appear in the coefficient Z in the way shown in (|34]) . 

In all cases relevant in the present context the matrix elements will be matrix elements 
involving hadronic states, which are in most cases impossible to calculate from first principles. 
However, eq.(^) allows to extract the short distance piece, i.e. the logarithms of the large mass 
rriq and to separate it into the Wilson coefficients. 

The anomalous dimension may be calculated in perturbation theory in powers of the cou- 
pling constant g of the theory. In general, in a renormalizable theory the coupling constant 
depends on the scale \x at which the theory is renormalized. The scale dependence of the 
coupling constant is determined by the f3 function 

g(j*)=PM- (38) 



rfln/i 

In a mass independent scheme the renormalization group functions 70 and (3 will depend on 
the scale \x only through their dependence on the coupling constant 

P = P(g(ji)) u = iMv)). (39) 

Hence we may rewrite the renormalization group equation fl5Bp as 




(^ + ^ + 7^)1^(^1=0. (40) 

The renormalization group functions (3 and 70 are calculated in perturbation theory; the first 
term of the j3 function on QCD is obtained from a one-loop calculation and is given by 

2 \ ,3 



P(9) = -j^y 2 {n--n f ) g * + ..., (41) 

where rif is the number of active flavors, i.e. the number of flavors with a mass less than tuq. 
With this input the renormalization group equation may be solved to yield 

where 71 is the first coefficient in the perturbative expansion of the anomalous dimension 
7j = li9 2 + ■ • • and a s (/x) is the one loop expression for the running coupling constant of QCD 

12tt 

as{fX) ~ (33-2n f )H^/Al CD ) (43) 




which is obtained from solving ( pq ) using (f4T|). This expression corresponds to a summation 
of the leading logarithms {a s \nmo) n which is achieved by a one-loop calculation of the renor- 
malization group functions (3 and 7q; in other words, in this way a resummation of the first 
column of the expansion (|3~4j) is obtained. 
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In a similar way one may also resum the second column of (^4|), if the renormalization group 
functions (3 and 7 are calculated to two loops and the finite terms of the one loop expression 
are included. 

Finally, the case we have considered as an example is indeed very simple; in general all op- 
erators of a given dimension may mix under renormalization, i.e. instead of a simple anomalous 
dimension a matrix of anomalous dimensions may occur. For more details on this I refer the 
reader to a textbook discussion of these issues as given e.g. in E| . 



2.3 Heavy Quark Symmetries 

The main impact of the heavy quark limit is due to two additional symmetries which are 
not present in full QCD; the first is a heavy flavour symmetry and the second one is the so- 
called spin symmetry. The presence of these symmetries implies Wigner-Eckart theorems for 
transition matrix elements which have far-reaching phenomenological consequences. 

We shall first study the heavy flavour symmetry. The interaction of the quarks with the 
gluons is flavour independent; all flavour dependence in QCD is only due to the different quark 
masses. In the 1 / uiq expansion the leading order Lagrangian is mass independent and hence a 
flavour symmetry appears relating heavy quarks moving with the same velocity. 

For the case of two heavy flavours b and c one has to leading order the Lagrangian || 

£ heavy = K(v ■ D)b v + C v (v ■ D)c v , (44) 

where b v (c v ) is the field operator h v for the b (c) quark moving with velocity v and D = d + igA 
is the QCD covariant derivative. This Lagrangian is obviously invariant under the SU(2)hf 
rotations 

I ( [' ) U v eSU(2) HF . (45) 

We have put a subscript v for the transformation matrix U, since this symmetry only relates 
heavy quarks moving with the same velocity. 

The second symmetry is the heavy-quark spin symmetry. As is clear form the Lagrangian 
in the heavy- mass limit, both spin degrees of freedom of the heavy quark couple in the same 
way to the gluons; we may rewrite the leading-order Lagrangian as 

£ = h+ s (ivD)h+ s + h- s (ivD)h-\ (46) 

where hf s are the projections of the heavy quark field on a definite spin direction s 

ht s = ~(l± l5 j)h v , s-v = 0. (47) 

This Lagrangian has a symmetry under the rotations of the heavy quark spin and hence all the 
heavy hadron states moving with the velocity v fall into spin-symmetry doublets as itlq — > oo. 
In Hilbert space this symmetry is generated by operators S v (e) as 

[h v ,S v (e)]=itf-y s h v (48) 

where e with e 2 = — 1 is the rotation axis. The simplest spin-symmetry doublet in the mesonic 
case consists of the pseudoscalar meson H(v) and the corresponding vector meson H*(v,e), 
since a spin rotation yields 

exp (<5„(e)|) \H(v)} = (-i)\H*(v, e)>, (49) 
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where we have chosen an arbitrary phase to be (— i). 

In the heavy-mass limit the spin symmetry partners have to be degenerate and their splitting 
has to scale as 1/toq. In other words, the quantity 

\ 2 = ~(M 2 H *-M 2 H ) (50) 

has to be the same for all spin symmetry doublets of heavy ground state mesons. This is 
well supported by data: For both the (£>,£>*) and the (D,D*) doublets one finds a value of 
A2 ~ 0.12 GeV 2 . This shows that the spin-symmetry partners become degenerate in the infinite 
mass limit and the splitting between them scales as l/m,Q. 

In the infinite mass limit the symmetries imply relations between matrix elements involving 
heavy quarks. For a transition between heavy ground-state mesons H (either pseudoscalar or 
vector) with heavy flavour / (/') moving with velocities v (v'), one obtains in the heavy-quark 
limit 

(H^\v')\h { pTh^\H^\v)) = £(W) Tr {H&)TH(v)} , (51) 

where T is some arbitrary Dirac matrix and H(v) are the representation matrices for the two 
possibilities of coupling the heavy quark spin to the spin of the light degrees of freedom, which 
are in a spin- 1/2 state for ground state mesons 



(1 + fyj 5 , (qQ) meson 
U(v) = { (1 + f)i I", (qQ) meson (52) 

with polarization e. 

Due to the spin and flavour independence of the heavy mass limit the Isgur-Wise function £ is 
the only non-perturbative information needed to describe all heavy to heavy transitions within 
a spin-flavour symmetry multiplet. 



Excited mesons have been studied in |25| . They may be classified by the angular momentum 
of the light degrees of freedom j t , which is coupled with the heavy quark spin S to the total 
angular momentum J of the meson. Furthermore, the orbital angular momentum £ determines 
the parity P = (— 1) £+1 of the meson. For a given I > we can have ji — £ ± 1/2 and the 
coupling of the heavy quark spin yields two spin symmetry doublets (J = £ — 1, J = £) and 
( J = £, J = £+ 1). For example, the lowest positive parity 1 = 1 mesons are two spin symmetry 
doublets (0 + , 1 + ) and (l + ,2 + ). In the D meson system these states have been observed [^(J 



and behave as predicted by heavy quark symmetry 

Similarly as for the mesons heavy-quark symmetries imply that only one form factor is 
needed to describe heavy to heavy transitions within a spin flavour symmetry multiplet; in 
other words, there is an Isgur Wise function for each multiplet. 

The ground state baryons have been studied in [|8|, |29|, |30| . According to the particle data 
group they are classified as follows 

Ah = [(qq')oh} 1/2 ~' h = [{qs) h] 1/2 (53) 
= [(<?</) 1^] 1/2 = [{qs)ih) 1/2 £L h = [(ss)i/i]i /2 (54) 

K = KwVls/a 2; = [(g S )!% 2 n; = [(ss)i% 2 . (55) 

Here, q, q' refer to u and d quarks, q ^ q' for the but q may be the same as q' for the 
and The first subscript (0, 1) is the total spin of the light degrees of freedom, while the 
second subscript (1/2, 3/2) is the total spin of the baryon. 
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Spin symmetry forces these baryons into spin symmetry doublets. For the A-type baryons 
(|53| ) the spin rotations are simply a subset of the Lorentz transformations, since the light degrees 
of freedom are in a spin-0 state. The corresponding spin symmetry doublet is in this case given 
by the two polarization directions of the heavy baryon. From the point of view of heavy quark 
symmetries the A-type baryons are the simplest hadrons, although from the quark model point 
of view they are composed of three quarks. 

The baryons with the light degrees of freedom in a spin one state may be represented by a 
pseudovector-spinor object i? M with v^R^ 1 = []. In general ■y^R^ 7^ because i? M contains spin 
1/2 contributions as well as spin 3/2 parts. In other words, i? M contains a Rarita-Schwinger 
field as well as a Dirac field. Under Lorentz transformations _R M behaves as 

R»{v) -> A^ u D(A)R u (Av), (56) 

where A^ u and D(A) are the Lorentz transformations in the vector and spinor representation 
respectively, while under spin rotations we have 

R"(v) - -^fiR^v). (57) 

The spin-3/2 component of the the pseudovector-spinor object corresponding to the Ej* is 
projected out by contracting with 7^ 

7^* = 0. (58) 
The rest of the independent components of R correspond to E/i baryon: 

^ = "^(7^ + ^)75^, (59) 

where u^ h is the Dirac spinor of the E^ state. Similar expressions hold for the non-strange 
baryons and f^*\ 



The spin rotation (57) transform the E-like baryons into the E* states and vice versa. Thus 
the spin symmetry doublets for the ground state baryons are given by the two polarization 
directions of the baryons in fl5"3"D , and by the two states with corresponding light quark flavour 
numbers in fl54f) and ([55|). 

Similar to the case of mesons one may derive a Wigner-Eckart theorem for the spin symmetry 
doublets of the baryons 

(A h (v)\hrti\A h ,(v')) = A(v ■ v')u ih {v)Tu ihl {v'), (60) 

where we have allowed for the possibility of two heavy quark flavours h and h'. In the same 
way, one obtains two form factors for the E^ — > E^. 

(^\v)\h v Th v/ \^) (61) 
= R% (v)TR u M (v) \B(v ■ v')g^ + C(v ■ v')v'v v . 

Finally, parity does not allow for transitions between A and £**) type baryons 

(p£\v)\h v Th v ,\A h (v)) = 0, (62) 

1 One could as well represent the light degrees of freedom by an antisymmetric tensor instead of a pseudovec- 
tor; this is a completely equivalent formulation p0[ ]. 
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and hence these transitions are not only suppressed by the flavour symmetry of the light degrees 
of freedom, but additionally by heavy quark symmetry. 

Excited baryons my be studied along the same lines as for the mesons. The spin symmetry 
doublets as well as the restrictions on transition matrix elements have been studied in J25 . 



Heavy quark symmetries thus lead to a strong reduction of the number of independent from 
factors that describe current induced transitions among heavy hadrons. In addition to that the 
symmetries even allow us to obtain the normalization of some of these form factors. Since the 
currents 

J hh> = h vl X = vji v hf v (63) 
are the generators of heavy flavour symmetry in the velocity sector v , the normalization of 



the Wigner-Eckard theorems (pi , 60 J6T|) is known at the non-recoil point v = v' . By standard 



arguments one obtains for the mesons 

CK = 1) = 1, (64) 

while the corresponding relation for the baryons is 

A(vv' = 1) = ^mK h mK h , (65) 
B(vv' = 1) = Jm (,)m (,), (66) 

h h' 



where the factor involving the square root of the masses means that the hadron states in (|60| ) 
are normalized relativistically. 

Up to now we have considered only the consequences of heavy quark symmetries for the 
leading terms of the I/tuq expansion. However, the additional symmetries also restrict the 
subleading terms and one of these restrictions is called Lukes theorem |TT| . It is a generalization 



of the Ademollo Gatto theorem fl32| , which states that in the presence of explicit symmetry 
breaking the matrix elements of the currents that generate the symmetry are still normalized 
up to terms which are second order in the symmetry breaking interaction. 

For the case at hand the relevant symmetry is the heavy flavor symmetry. This symmetry 
is an SU (2) symmetry and is generated by three operators Q± and Q3 with 

Q+ = J d 3 xb v (x)c v (x) Q- = J d 3 x c v (x)b v (x) 

Qi = J d 3 x(b v (x)b v (x) - c v (x)c v (x)) (67) 
[Q+,Q-] = Q* [Q + ,Qs} = -2Q + (Q+V = Q- 

Let us denote the ground state flavour symmetry multiplet as \B) and \D). Then the 
operators act in the following way 

Q 3 \B) = \B) Q 3 \D) = -\D) 

Q+\D) = \B) Q_\B) = \D) (68) 
Q+\B) = Q-\D) = 0. 

The Hamiltonian of this system has a I/tuq expansion of the form 



H = + HP + -^H[ b} + -^H[ c) + ■■■ (69) 
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Hgyrarn ^break- 



In the second equation, the first line is still symmetric under heavy flavour 577(2) while the 
term in the second line does not commute any more with Q±, but it still commutes with Q%. 
In other words, to order I/tuq we still have common eigenstates of H and Q3, which we shall 
denote as \B) and \D). Sandwiching the commutation relation we get 

1 = (B\Q 3 \B) = (B\[Q + ,Q_]\B) (70) 
= [(B\Q+\n)(n\Q-\B) - (B\Q_\n)(n\Q + \B) 



^[\(B\Q + \n)\ a -\(B\Q.\n)\ 



where \n) form a complete set of states of the Hamiltonian Hsymm+H^^- The matrix elements 
may be written as 

(B\Q±\n) = 1 (B\[H break ,Q ± ]\n) (71) 

where Eb and E n are the energies of the states \B) and |n) respectively. In the case \n) = \D) 
the matrix element will be of order unity, since both the numerator as well as the energy 
difference in the denominator are of the order of the symmetry breaking. For all other states 
the energy difference in the denominator is non- vanishing in the symmetry limit, and hence 
this difference is of order unity; thus the matrix element for these states will be of the order of 
the symmetry breaking. From this we conclude 



(B\Q+\D) = 1 + 



1 



2rrih 2m 



(72) 



In particular, the weak transition currents at the non-recoil point v = v ' are proportional to 
these symmetry generators and hence we may conclude that for some of these matrix elements 
we only have corrections of the order 1/ttiq. 

Another restriction on the 1/rriQ expansion is imposed by the so-called reparametrization 
invar iance which is basically the remnant of the original Lorentz covariance of full QCD. 
The full theory depends only on the momentum of the heavy quark P, and the splitting of 
this momentum in an on-shell part w,qV and a residual momentum k corresponding to the 
covariant derivative acting on the heavy static field h v is arbitrary. Formally this means that 
the Lagrangian does not depend on the velocity v, if all orders of the 1/toq expansion are 
included; the v dependence only enters once the expansion is truncated. 

Using the representation flT9| ) and (|IBD the Lagrangian is invariant under the transformation 

v — > v + Sv v ■ 8v = (73) 

h v ^h v + ^(l + P- 1 -— F ,tp) K 

2 \ 2m + iv D ) 

iD — > —m5v. 

This invariance is the so-called reparametrization invariance, which has non-trivial conse- 
quences, since it relates terms of different orders of the 1/m expansion. 
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3 Application to Exclusive Decays 



The heavy mass limit and the resulting additional symmetries allow us to restrict the matrix 
elements which occur in weak transitions of heavy hadrons. We shall consider in the following 
in some detail the semileptonic b — > c transition, which we shall treat as a heavy to heavy 
decay. In section 3.2 we investigate the consequences of heavy quark symmetries for transitions 
of the heavy to light type. 



3.1 Transitions of the type Heavy — ► Heavy 



For the case of a heavy to heavy transition the Wigner Eckart theorem (|5l|) implies that there 
is only a single form factor which describe the weak decays of heavy hadrons; furthermore, the 
heavy mass limit yields the normalization of this form factor at the kinematic point v = v'. 

Treating both the b and the c quark as heavy, the semileptonic decays B — > D^*'£u are the 
phenomenologically relevant examples. The matrix elements for these transitions are in general 
parametrized in terms of six form factors 



(D(v')\c^b\B(v)) = + <) + £-(y)(^ - <) 

(D m (v',e)\cry,J>\B(v)) = i^/m B m D »i v {y)e m p P e* a v' p v p 
(D*(v',e)\cj^ 5 b\B(v)) = Vm fl m fl . [Ui(v)(w' + l)e* - £ M (y)(e*v)v 

-6i2(2/)(^)< 



(74) 
(75) 

(76) 



where we have defined y = vv'. Due to the Wigner Eckart theorem (|5l|) these six from factors 
are related to the Isgur Wise function by 



ay) = Z(y) for i = +,v, Ai, A3, &(y) = o for i = -, A2. 



(77) 



Since heavy quark symmetries also yield the normalization of the Isgur Wise function, we know 
the absolute value of the differential rate at the point v — v' in terms of the meson masses 
and V c b. Hence we may use this to extract V c b from these decays in a model independent 
way by extrapolating the lepton spectrum to the kinematic endpoint v = v'. Using the mode 
B — > D^tv one obtains the relation 



lim 



1 



dY 



G 



sj{vv') 2 - 1 d{vv') 4vr 3 



^V ch \ 2 {m B - m D *) 2 m 3 Dt \£ A1 (l) 



(78) 



where £ai is equal to the Isgur Wise function in the heavy mass limit, and hence £ai(1) = 1- 

Corrections to this relation have been calculated along the lines outlined above in leading 
and subleading order. A complete discussion may be found in more extensive review articles 
(see e.g. Neubert's review 0), including reference to the original papers. Here we only state 
the final result 



Wl)=* 6 / 25 
25 14 



+z < x 

\54 27 



1 + 1.561 

9/25 



a s {m c ) - a s {m b ) 8a s (m c ) 
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3tt 



(79) 
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where we use the abbreviations 



a,\m r 



x 



m b 



a s {m b y 

and m is a scale somewhere between m b and m c . 

Up to the term 5i/ m 2 all these contributions may be calculated perturbatively, including 
the dependence on z. The quantity b\i m 2 parametrizes the non-perturbative contributions, 
which enter here at order 1/m 2 . These corrections may be expressed in terms of the kinetic 
energy Ai, the chromomagnetic moment A2, which are given in 



and (28) respectively, and 
matrix elements involving time-ordered products between the current and the corrections of 
the Lagrangian 



/m 2 



/ 1 



\2m c 



1 



-Ai + A 2 



(80) 



d 4 xd*y (B*(v,e)\T \C { b 1] {x)b vCv C^ (y)] \D*(v,e)) 



v ' 2^M B M D 
+ 0(l/ml,l/m 2 b ,l/(m c m b )), 

where Cq* is the first order Lagangian for the quark Q as given in ([H]) or ( p3|) and Mb (Mp) 
are the masses of the B (D) meson in the heavy quark limit. Here we display only the largest 
contribution of order the complete expression, including the 1/m 2 and l/(m c m b ) terms, 

may be found in 

Thus the correction S 1 / m 2 is given in terms of Ai defined in (0), A 2 given in (|28"D and 
a non-local matrix element involving a time-ordered product. The problem concerning the 
determination of Ai has been considered already above; similarly it is not easy to obtain infor- 
mation on the matrix element involving the time-ordered product, and thus the corrections of 
order 1/m 2 will finally limit our ability to determine the CKM matrix element V cb in a model 
independent way, at least using the approach described here. 

Various estimates for 5i/ m 2 have been given in the literature. The first estimate of this 
correction has been given in |}4] using the GISW model [Q, which is based on a wave function 
for the light quark. In this work 5 m 2 = —2% ... — 3% has been obtained. Another estimate 

but both estimates have been criticized 



35 



with weaker assumptions yields S m 2 = ... — 5% 
recently as being too small. Based on heavy flavour sum rules it has been argued in [37 that the 
1/m 2 corrections can be quite large 5 m 2 = 0% ... — 8% [|3"?|| . These various estimates indicate 
the size of the theoretical error involved in the determination of V cb from the exclusive channel 
B — > D*£uf, a generally accepted value for these corrections has been given recently p8 



-0.055 ±0.025 



(81) 



from which one obtains 



=0.91 ±0.03 (82) 

This result has been used to extract V cb from CLEO |39| as well as from LEP data |I(J. The 
values obtained are 



\V cb \ 

\Vcb\ 



0.0386 ±0.0019 ±0.0020 ±0.0014 CLEO 
0.0392 ± 0.0025 ± 0.0027 ± 0.0015 ALEPH 



(83) 
(84) 



where ([82]) has been used. Note that the third error in \ V cb \ is due to the theoretical uncertain- 
ties, which by now almost match the experimental ones. 
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3.2 Transitions of the type Heavy — > Light 

Heavy quark symmetries may also be used to restrict the independent form factors appearing 
in heavy to light decays. For the decays of heavy mesons into light 0~ and 1~ particles heavy 
quark symmetries restrict the number of independent form factors to six, which is just the 
number needed to parametrize the semileptonic decays of this type. Furthermore, no absolute 
normalization of form factors may be obtained from heavy quark symmetries in the heavy to 
light case; only the relative normalization of B meson decays heavy to light transitions may be 
obtained from the corresponding D decays. 

In general we shall discuss matrix elements of a heavy to light current which have the 
following structure 

J = (A\iVh v \H(v)), (85) 

where T is an arbitrary Dirac matrix, I is a light quark (u, d or s) and A is a state involving 
only light degrees of freedom. 

Spin symmetry implies that the heavy quark index hooks directly the to the heavy quark 
index of the Dirac matrix of the current. Thus one may write for the transition matrix ele- 
ment (0) 

(A\£Th v \H{v)) = Tr {M A TH{v)) (86) 

where the matrix H(v) representing the heavy meson has been given in (^). The matrix M.a 
describes the light degrees of freedom and is the most general matrix which may be formed from 
the kinematical variables involved. Furthermore, if the energies of the particles in the state A 
are small, i.e. of the order of Aqcd, the matrix M.a does not depend on the heavy quark; in 
particular it does not depend on the heavy mass ran- In the following we shall discuss some 
examples. 

The first example is the heavy meson decay constant, where the state A is simply the 
vacuum state. The heavy meson decay constant is defined by 

(0\h^ 5 h v \H{v)) = f H m H v^ (87) 

and since \A) = |0) the matrix M.® is simply the unit matrix times a dimensionful constant^ 
and one has, using ( |86D 

Wl^h v \H{v)) = ft Tr (jj b H(v)) = 2«v^y ( 88 ) 

As discussed above the constant k does not depend on the heavy mass and thus one infers the 
well-known scaling law for the heavy meson decay constant from the last two equation 

fn oc — L (89) 

Including the leading and subleading QCD radiative corrections one obtains a relation between 
f B and f D 



m b \a s (m c ) 



71 



f D ~ 0.69/d. (90) 



2 Note that contributions proportional to f may be eliminated using 

H{v)i> = -H(v). 
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The second example are transitions of a heavy meson into a light pseudoscalar meson, which 
we shall denote as tt. The matrix element corresponding to (|55D is 

J P =(7T(p)\iTh v \H(v)}, (91) 

where p is the momentum of the light quark, 

The Dirac marix Aip for the light degrees of freedom appearing now in (|86|) depends on 
p and v . It may be expanded in terms of the sixteen independent Dirac matrices 1, 75, 7^, 
757^, and a^ v taking into account that it has to behave like a pseudoscalar. The form factors 
appearing in the decomposition of Aip depend on the variable v ■ p, the energy of the light 
meson in the rest frame of the heavy one. In order to compare different heavy to light transition 
by employing heavy flavor symmetry this energy must be sufficiently small, since the typical 
scale for the light degrees of freedom has to be of the order of Aqcd to apply heavy quark 
symmetry^. For the case of a light pseudoscalar meson the most general decomposition of Aip 
is 

Aip = ^TpA{r,) l5 + ——B{ri)7sh (92) 
where we have defined the dimensionless variable 

(93) 

A QCD 

The form factors A and B are universal in the kinematic range of small energy of the 
light meson, i.e. where the momentum transfer to the light degrees of freedom is of the order 
Aqcd] in this region rj is of order unity. This universality of the form factors may be used to 
relate various kinds of heavy to light transitions, e.g. the semileptonic decays like D — > Treu, 
D — > Kev or B — > -Rev and also the rare decays like B — > K£ + £~ or B — > tt£ + £~ where £ 
denotes an electron or a muon. 

As an example we give the relations between exclusive semileptonic heavy to light decays. 
The relvant hadronic current for this case may be expressed in terms of two form factors 

(7r(p)|l 7 (l - *y s )hv\H(v)) = F x (v ■ p)m H v^ + F 2 (v ■ p) Pfl (94) 

= F + (vp){m H v tl +p tl ) + F_(vp)q il 



where 



F ± {vp)= l -{F l {vp)±F 2 {vp)) (95) 



2 

Inserting this into (|9TD one may express F± in terms of the universal form factors A and B 



FAv-p) = F+(v ■ p) + FJv ■ p) =-2i/— Afa) (96) 

V m H 



F 2 {v ■ p) = F + {v ■p)-F4vp) = -2J^B(r]) (97) 
From these relations one may read off the scaling of the form factors with the heavy mass which 



was already derived in |4l 



3 Note that in this case the variable v -p ranges between and m#/2 where we have neglected the pion mass. 
Thus at the upper end of phase space the variable v ■ p scales with the heavy mass and heavy quark symmetries 
are not applicable any more. 
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This may be used to normalize the semileptonic B decays into light mesons relative to the 
semileptonic D decays. One obtains 



Note that F + for the B decay is expressed in terms of F + and F_ for the D decays. In the 
limit of vanishing fermion masses only F + contributes, which means that the _F_ contribution 
to the rate is of the order of mi ept on/mH- Thus it will be extremely difficult to determine 
experimentally. 

The case of a heavy meson decaying into a light vector meson may be treated similarly. 
The matrix element for the transition of a heavy meson into a light vector meson (denoted 
generically as p in the following) is given again by (|85f) and is in this case 

J v = (p(p,e)\£Th v \H(v)). (99) 

Using fl86|) one has 

(p(p,e)\£Th v \H{v)) = Tr (M v TH(v)), (100) 

where now the Dirac matrix M.y has to be a linear function of the polarization of the light 
vector meson. 

The most general decomposition is given in terms of four dimensionless form factors 

M v = ^TpC(r))(v ■ e) + — == D{rj){v ■ e){> + Jv^E{-q)<f. + -^=F( V )^ (101) 

where the variable r\ has been defined in (|93|). 

Similar to the case of the decays into a light pseudoscalar meson ( |100| ) may be used to 



relate various exclusive heavy to light processes in the kinematic range where the energy of the 
outgoing vector meson is small. For example, the semileptonic decays D — ► pev, D — ► K*ev 
and B — > pev are related among themselves and all of them may be related to the rare heavy 
to light decays B -> K*£ + £~ and B -> p£ + £~ with £ = e, p. 

Data on these decays are still very sparse; there are first measurements of the decays B — > 
k£v and B — ► p£t/ from CLEO P2| |, from which total rates may be obtained. From this one 



may extract a value of V u b by employing form factor models, and the value given by CLEO is 

\V ub \ = (3.3 ± 0.2+°;3 ± 0.7) x 10~ 3 (102) 

where the last uncertainty represents the variation of the result between different models. In 
order to perform a model independent determination along the lines discussed above a good 
measurement of the lepton energy spectra in these decays is needed. 

Finally we comment on the heavy to light transitions of baryons. For the A-type heavy 
baryons Q55p spin symmetry relates different polarizations of the same particle and thus imposes 
interesting constraints. Consider for example the matrix element of an operator £Th v between 
a heavy Aq and a light spin-1/2 baryon Bi. It is described by only two form factors, 

(Bt(p)\iTh v \K Q {v)) =u l (p){F 1 {vp)+^F 2 {vp)}Tu AQ {v). (103) 

Thus in this particular case spin symmetry greatly reduces the number of independent Lorentz- 
invariant amplitudes which describe the heavy to light transitions. 
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This has some interesting implications for exclusive semileptonic A c decays. For the case of a 
left handed current T = 7^(1 — 75), the semileptonic decay A c — > A£p£ is in general parametrized 
in terms of six form factors 

(A(p)|$7„(l - 1b)c\A c (v)) = u(p) [/ l7/i + ih<J»»q v + / 3 g1 u{p') 

+ u(p) [#17/, + m^^q v + g-^} j5u(p'), (104) 

where p' = m Ac v is the momentum of the A c whereas q = m\ c v — p is the momentum transfer. 
From this one defines the ratio Ga/Gv by 

Ga _ 9^ = 0) 

G v - fl ( q i = oy (105) 

In the heavy c quark limit one may relate the six form factors fi and gi (i — 1, 2, 3) to the 
two form factors Fj (j = 1,2) 

fi = -9i = F, + ^F 2 (106) 

/ 2 = h = -92 = -93 = —F 2 (107) 

from which one reads off Ga/Gv = — 1- This ratio is acessible by measuring in semileptonic 
decays A c — > \lvi the polarization variable a 

2G A G V 

which is predicted to be a = —1 in the heavy c quark limit. The subleading corrections to the 



heavy c quark limit have been estimated and found to be small [H3 



a < -0.95, (109) 

and recent measurements yield 

a = -0.91 ± 0.49 ARGUS [44] (110) 

a = -0.89^;™ CLEOH (111) 



and are in satisfactory agreement with the theoretical predictions. 

Recently the CLEO collaboration also measured the ratio of the form factors iq and _F 2 , 
averaged over phase space. Heavy quark symmetries do not fix this form factor ratio, at least 
not for a heavy to light decay, while for a heavy to heavy decay the form factor F 2 vanishes in 
the heavy mass limit for the final state quark. CLEO measures [f|(J 

^\ = -0.25 ±0.14 ±0.08 (112) 

F 2 1 phase space 



which is in good agreement with model estimates [4/] 
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4 The 1/rriQ Expansion in Inclusive Decays 



For inclusive decays a I/tuq expansion is obtained for the rates by an approach similar to 



the one known from deep inelastic scattering [p||- ||13| . The first step consists of an operator 
product expansion (OPE) which yields an infinite sum of operators with increasing dimension. 
The dimensions of the operators are compensated by inverse powers of a large scale, which 
is in general of the order of the heavy mass scale. The decay probability is then given as 
forward matrix elements of these operators between the state of the decaying heavy hadron; 
these matrix elements still have a mass dependence, which then may be extracted in terms of 
a 1/ rriQ expansion using HQET as for exclusive decays. 

The method described below also allows us to deal with inclusive non-leptonic processes and 
hence in principle opens the possibility for a calculation of lifetimes and branching fractions 
in the framework of the I/tuq expansion. This is remarkable, since non-leptonic processes are 
usually very hard to deal with, in particular the 1/m.Q expansion has not (yet ?) brought any 
success in the field of exclusive non-leptonic decays. 

Applying the OPE to the energy spectra of the charged lepton in inclusive semileptonic 
decays of heavy mesons, the relevant expansion parameter is not 1/mq, but rather l/(mg— 2Eg); 
the denominator is thus the energy release of the decay. In almost all phase space the energy 
release is of the order of the heavy mass; it is only in the endpoint region that it becomes small 
and hence the expansion breaks down. This problem may be fixed by a resummation of terms 
in the operator product expansion, which strongly resembles the summation corresponding to 
leading twist in deep inelastic scattering. Analogously to the parton-distribution function, a 
universal function appears, which determines all inclusive heavy-to-light decays. 

4.1 Operator Product Expansion for Inclusive Decays 

The effective Hamiltonian for a decay of a heavy (down-type) quark is in general linear in the 
decaying heavy flavoured quark 

H eff = QR (113) 

where the operator R describes the decay products. In the following we shall consider semilep- 
tonic decays, for which 

R sl = ^V Qq 7 „(1 - 75 )g (^7 M (1 - 75)*), (114) 

where q is an up-type quark (c or u, since we shall consider b decays). Similarly, for non-leptonic 
decays the Cabbibo allowed contribution corresponds to 

Rm = §7^V Qq V* q „ [(C+(m b ) + C_(m 6 )) 7M (l - l,)q (gVU - Ts)?") (115) 

+ (C + (m b ) - C_(m 6 )h M (l - IsW (gV(l - 7fe)ff)] , 

where q' (q") is a down-type (up-type) quark and Vq q the corresponding CKM matrix element. 
The coefficients C±(?7i&) are the QCD corrections obtained from the renormalization group 



running between M\y and mj; in leading logarithmic approximation these coefficients are [48 



C±(m b ) = [^t^p, with 7+ = 6 - = -V (116) 

V ' [ a s {m 2 b ) ] ' I+ 33-2N f 2 1 K ' 
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where a s (/i) is the onle-loop expression fl43[) for the running coupling coupling constant of QCD. 
Finally, for radiative rare decays we have 

R rare = -^V tb V*C 7 (m b )—m b a^(l + l5 )sF^ (117) 

where Cr(rrib) is again a coefficient obtained from running between Mw and rrib. Its value is 
CV(m&) ~ 0.3, the corresponding analytical expression may be found in PS[| . 

The inclusive decay rate for a heavy hadron H containing the quark Q may be related to a 
forward matrix element by 

T^Y,(^ 4 (PB-Px)\(X\H eff \H(v)}\ 2 = [ d*x(H(v)\H eff (x)Ht ff (0)\H(v)) (118) 
x J 

= 2 Im / d*x(H(v)\T{H ef f(x)Hl f (0)}\H(v)). 



where \X) is the final state, which is summed over to obtain the inclusive rate. 

The matrix element appearing in ( |118| ) contains a large scale, namely the mass of the heavy 
quark. The first step towards a 1/rriQ expansion is to make this large scale explicit. This may 
be done by a phase redefinition. This leads to 

r oc 2 Im | d A xe- im Q vx (H(v)\T{H eff (x)Hl ff (Q)}\H(v)) (119) 

where 

H ef f = Q V R Q v = e- m ^ vx Q (120) 

This relation exhibits the similarity between the cross-section calculation in deep inelastic 
scattering and the present approach to total rates. In deep inelastic scattering there appears a 
large scale, which is the momentum transfer to the leptons, while here the mass of the heavy 
quark appears as a large scale. 

The next step is to perform an operator product expansion of the product of the two 
Hamiltonians. After the phase redefinition the remaining matrix element does not involve large 
momenta of the order of the heavy quark mass any more and hence a short- distance expansion 
becomes useful, if the mass mg is large compared to the scale A determining the matrix element. 
The next step is thus to perform an operator-product expansion, which has the general form 



J d * xe *m QVX {H (v)\T{n eff (x)Hl ff (0)}\H(v)} = £ ( — j C n+ M{H{v)\O n+ ,\H{v))„ 



where O n are operators of dimension n, with their matrix elements renormalized at scale /i, 
and C n are the corresponding Wilson coefficients. These coefficients encode the short distance 
physics related to the heavy quark mass scale and may be calculated in perturbation theory. 
All long distance contributions connected to the hadronic scale A are contained in the matrix 
elements of the operators O n+ z- 

Still the matrix elements of O n+ s are not independent of the heavy quark mass scale, but 
this mass dependence may be expanded in powers of 1/rriQ by means of heavy quark effective 
theory. This is achived by expanding the heavy quark fields appearing in the operators O n using 



(18) (or, equivalently, (|22|) ) as well as the states by including the corrections to the Lagrangian 
given in (|19D (or (^3|)) as time-ordered products. In this way the mass dependence of the total 
decay rate may be accessed completely within an expansion in I/tuq. 
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The lowest-order term of the operator product expansion are the dimension-3 operators. 
Due to Lorentz invariance and parity there are only two combinations which may appear, 
namely Q v fQ v or Q V Q V . Note that the Q v operators differ from the full QCD operators only 
by a phase redefinition, and hence Q V $Q V = Q^Q and Q V Q V = QQ. The first combination is 
proportional to the Q-number current Qj^Q, which is normalized even in full QCD, while the 
second differs from the first one only by terms of order 1/toq 



(iDf 



[ivDY + 



-a llu G tlv 



h v + 0{l/m%) 



(121) 



where G pu is the gluon field strength. 

Thus the matrix elements of the dimension-3 contribution is known to be normalized; in 
the standard normalization of the states this implies 



>H(v)\0 3 \H(v)) = {H{v)\Q v fH v \B{v)) = 2m 



H 



(122) 



where m# is the mass of the heavy hadron. To lowest order in the heavy mass expansion we 
may furthermore replace mg = tuq and hence we may evaluate the leading term in the l/m q 
expansion without any hadronic uncertainty. Generically the dimension-3 contribution yields 
the free quark decay rate. This has been previously used as a model for inclusive decays, but 
now it turns out to be the first term in a systematic l/mq expansion of total rates. 

A dimension-four operators contains an additional covariant derivative, and thus one has 
matrix elements of the type 



(H(v)\0,\H(v)) cc (H{v)\Q v TD fl Q v \H{v)) = A r v p 



(123) 



Since the equations of motion apply for this tree level matrix element, one finds that the 
constant Ar has to vanish, and thus there are no dimension-four contributions. This statement 
is completely equaivalent to Lukes theorem ||31|| , since we are considering a forward matrix 
element, i.e. a matrix element at zero recoil 
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The first non-trivial non-perturbative contribution comes from dimension-5 operators and 
are of order I/itlq. For mesonic decays there are only the two parameters Ax and A2 given in 
(|27|) and fl28), which correspond to matrix elements of the subleading terms of the Lagrangian. 
They parametrize the non-perturbative input in the order l/nig. For Ag-type baryons the 
parameter A2 vanishes due to heavy quark spin symmetry, while the kinetic energy parameter 
Ai is non-zero as well. In the framework of the 1/mg expansion this leads to a difference in 
lifetimes between mesons and baryons. 



4.2 Calculation of Total Decay Rates 

In this subsection we shall collect the results for the total rates including the first non-trivial 
non-perturbative correction. 

Inserting R s i as given in (|114 ) one obtains for the total inclusive semileptonic decay rate 
B -> XJLv 



G 2 F ml 

1927T 3 



V cb \ 2 



Ai \ ( m c 
2m 2 I Vmfe 



9A 2 , / m e 
2m 2 Wi, 



where the two /,• are phase-space functions 



fi(x) 



1 - 8x 2 + 8x b 



x 



24x 4 log x, 
5 



(124) 



(125) 



1 - -x 2 - 8x A + 8x 6 + -x 8 + 8x 4 log x. 
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The result for B — > X u tvi is obtained from (|124|) as the limit m c 

V c b — > V u b 



T(B XJu) 



G 2 F ml 

1927T 3 



Ai - 9A 2 
2mg 



and the replacement 



(126) 



As it has been discussed above, the leading non-perturbative corrections in (|124 ) and ( |1 26| ) are 
parametrized by Ai and A 2 . Estimates for these parameters have been discussed in section 2; 
in order to estimate the total effect of the non-perturbative effects we shall insert a range of 
values —0.3 > Ai > —0.6 GeV 2 ; from this we obtain 



Ai - 9A 2 
2m 2 b 



-(3---4)% 



(127) 



This means that the non-perturbative contributions are small, in particular compared to the 
perturbative ones, which have been calculated some time ago 0, ^1] . For the decay B — > X u £u e 
the lowest order QCD corrections are given by 



T(B -> XJP t ) 



G 2 F ml 

1927T 3 



V, 



ub 



+ 3tt VT 



7T 



o.85iKi,rr 



128) 



and thus the typical size of QCD radiative corrections is of the order of ten to twenty percent. 
Similarly, one obtains the result for non-leptonic decays as 



T(B - X c ) = 3 



G\m\ 

1927T 3 



\m b / 



1 + 



1 



2m\ 



(Ai - 9A 2 ) 



48A 2 / 3 



m c \ 1 



nib J 2m^ 



2 ^2 



(129) 



where the coefficients Ai are given by combination of the Wilson coefficents C±{nib) ( 116 ) 



A 1 = -[C 2 _(mb) + 2C 2 + (m b )}, 



A, 



6 



[C 2 Jm b ) - C 2 _(m b )}, 



(130) 



and f 3 (x) 



x 



2\3 



is another phase space function. Again the non-perturbative corrections 



turn out to be small, in the region of a few percent compared to the leading term, and the 
perturbative corrections turn out to be much larger than this. 

Finally, for the rare decay B — > X s 7 one may as well calculate the non-perturbative contri- 
bution in terms of Ai and A 2 . One obtains 



Y{B - X sl ) 



aGj 
16V 



mi\V ts V t * d \ 2 \C 7 (m b )\< 



1 + 



2m\ 



(Ai-9A 2 ) 



(131) 



and the relative size of the non-perturbative corrections is the same as in the B X u li>i 
decays. 

Typically the non-pertubative corrections are much smaller than the radiative corrections. 
The only exception is the endpoint region of lepton energy spectra which receives both large 
perturbative as well as non-perturbative corrections. However, this is only a small region in 
phase space and the corections to the total rates remain moderate. 



4.3 Lifetimes of Heavy Hadrons 

The subject of heavy hadron lifetimes is strongly related to non-leptonic processes, which have 
been considered already some time ago ||52||-||55||; however the application of the 1/niQ expansion 
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FIG 1 



Figure 1: Diagrams for non-leptonic decays of B mesons. The upper diagram corresponds to 
the leading term in the 1/m.g expansion, the middle one to weak anihilation, and the lower one 
to Pauli interference. Diagrams taken from pH . 



has turned many assumptions into quantitative arguments. As outlined in the last section the 
1/rriQ expansion allows to calculate total rates, even for non-leptonic processes, and hence a 
QCD-based calculation of lifetimes becomes possible. A recent review of this subject is given 
in g. 

Studying the formulae obtained in the 1/tuq expansion up to order 1/mq one finds that 
lifetime differences between B mesons do not occur up to this level of the expansion; in other 
words, any difference between the B and the B° or the B s lifetimes are induced by effects of 
the order l/m b . These effects are due to dimension six operators of the four-quark form 

6 = (bTq)(qTb). (132) 

In fig.0 the relevant diagramms are shown. The two lower ones yield dimension six operators 
of the from ( |132j ); comparing to the terminology used in phenomenological models the middle 
diagram is called weak annihilation (WA) and the lower one Pauli interference (PI). 

The four quark operators appearing at order l/m^ are usually estimated applying the 
vacuum insertion assumption, althought this procedure has been criticized recently [j57l| . It 
amounts to replace 

(B\(bTq)(qTb)\B) -> (B\(bTq)\0) (0\(qTb)\B) 

and hence the parameter entering the estimates of the lifetime differences is the decay constant 
f B - 



The WA piece has been considered in |53j po| and has been found to be small, of the order 
of one percent. The dominant contribution comes from the PI diagram. For the lifetime of the 



B one obtains ||56 



T(B- 



T l/m2 (B) + AT PI (B- 



(133) 



where the PI contribution is 



AT PI (B- 



G 2 F ml 

1927T 3 



V cb \ 2 ■ 24tt 



I 2 
M 2 



B 



[C 2 + (m b ) - C 2 (m b ) + —(C 2 Jm b ) + C 2 _{m b ))} (134) 



B 



N, 



c 



where C± have been given in (|116|) . QCD radiative corrections change the sign of the PI 
contribution; the constructive interference at the scale M\y is turned into a destructive one at 
the scale m b prolonging the lifetime of the B~ relative to the Bq. Using HQET the running 
of the coefficients below m b has also been calculated; it has been found that the destructive 
interference is amplified by the running below m b , since 



AT PI (B- 



G 2 F m\ 2 
192tt3 IM 



24vr 2 Jb 



M 2 B 



k- 4 [(C 2 + (m b )-C 2 _(m b ))K 9 / 2 , 

- 1 -{k'I 2 - l)(C 2 + (m b ) - C 2 (m b ))} 



C 2 Jm b )+C 2 (m b ) 



(135) 
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where 

k = [ as (f h ?} ] 1 ' b , with 6=11- -n F (136) 
a s {m 2 b ) 3 

and Hhad is some small hadronic scale, which may be defined e.g. by a s (fihad) = 1- Putting in 
numbers the result is 



T 



1 + 0.05-- J -£ — - (137) 

r(B d ) (200 MeV) 2 V ; 

Hence the 1/toq expansion (in combination with vacuum insertion) predicts a slightly longer 
lifetime for the charged B meson compared to the neutral one. However, this is based on the 
vacuum insertion assumption which has been reconsidered recently in ||57|| ; it is claimed that 



possible non-factorizable contributions may be enhanced by large prefactors thereby invalidating 

(EH. 

However, the estimate (|137|) is compatible with data; the latest compilation yields |59| 



-(B + ) 

1.019 ±0.048 (138) 



Applying the 1/toq expansion for the lifetimes also to the D system one finds large cor- 
rections. From the experimental side it is knwon that the lifetime differences are very large 
and that the naive parton model expectation fails by a large margin. This may be taken as an 
indication that the c quark mass is indeed not large enough to justify a I/tuq expansion, at 
least for the non-leptonic processes. 

The lifetimes of the two neutral B mesons are equal up to terms of order l/mf, hence one 
expects 

f{B d ) = f{B s ) (139) 

to a good accuracy. Here f denotes the average lifetime of the two mass eigenstates of the 
B° — B° system. While the lifetimes in the B d system are practically the same such that the 
difference may be neglected, it has been pointed out that this is not necessarily the case in the 



B s system; the lifetime difference has been estimated in [60] to be 



AT(B S ) _ T(B s>short ) —T(B St i ong ) _ f j 



0.18 •- ^ — - (140) 

F(B S ) T(B S ) (200 MeV) 2 V ; 

Hence this difference may be as large twenty percent and thus it cannot be neglected any more 
e.g. in an analysis of B s -B s mixing. 

Finally we shall consider the b baryon lifetimes. Here we expect the differences between the 
meson and the baryon lifetimes to be of order 1/m 2 , since the matrix elements of the kinetic en- 
ergy operator as well as of the chromomagnetic moment operator are different between baryons 
and mesons; in particular, due to spin symmetry the matrix element of the chromomagnetic 
moment operator vanishes for A-type baryons, since the light degrees of freedom are in a spin-0 
state. Probably only the A-type baryons will be stable against electromagnetic and strong 
decays, and hence their lifetime will be determined by the weak decay of the heavy quark. 

Theoretical estimates of the A^ lifetime have been attempted in J56 , |57| and yield the 
expectation that 

r(A 6 ) > 



i Bi) 0A (141) 
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This has to be confronted with the recent data 



59 



= 0.80 ± 0.05. (142) 

Thus the lifetime is slightly below the expectation, although this is not yet a significant deviation 
from the prediction. In particular, older data indicated that this lifetime ratio could have been 
as low as 0.7; such a low value would clearly indicate a theoretical problem in the 1/toq 
expansion of the lifetimes, but the new data are in better agreement with the theoretical 
expectations. 



4.4 Lepton Energy Spectra 

The method of the operator-product expansion may also be used to obtain the non-perturbative 
corrections to the charged lepton energy spectrum f§ . In this case the operator product expan- 
sion is applied not to the full effective Hamiltonian, but rather only to the hadronic currents. 
The rate is written as a product of the hadronic and leptonic tensor 



dT 



Gl ' l V Qq \ 2 W, u A" u d(PS), 



Am 



(143) 



B 



where d(PS) is the phase-space differential. The short-distance expansion is then performed 
for the two currents appearing in the hadronic tensor. Redefining the phase of the heavy-quark 
fields as in (|ll|) one finds that the momentum transfer variable relevant for the short-distance 
expansion is txiqV — q, where q is the momentum transfer to the leptons. 

The structure of the expansion for the spectrum is identical to the one of the total rate. The 
contribution of the dimension-3 operators yields the free-quark decay spectrum, there are no 
contributions from dimension-4 operators, and the 1 / ml corrections are parametrized in terms 
of Ai and A 2 . Calculating the spectrum for B — > X c tv yields 



dT 

dy 



G 2 F \V cb \ 2 


ml 


1927T 3 




+ A ' 




[m b (l - 


v)\ 


3A 2 




m 2 b (l - 





+^[(6 + 5y) 



9(1 -y 



AR 



p)y 2 [{3(l-p)(l 
Ai 



R 2 ) - 2y(l - R 3 )} 



(144) 



y) 



(R 2 - 2R A 



Ai 



+ -^ 1 [5y-2(3-p)R 2 + AR- i 



3ml 



mi 



12R - (9 - 5p)R 2 + 10R 3 



+ o 



(A/[m b (l-y)}f 



where we have defined 



and 



m c \- 
mj 



R 



y 



y = 2E e / m b 

is the rescaled energy of the charged lepton. 

This expression is somewhat complicated, but it simplifies for the decay B 
then the mass of the quark in the final state may be neglected. One finds 



(145) 



(146) 



X„£i> since 



dT 

dy 



r<2 



Vub\ 2 m 

1927T 3 



2y 2 {2, - 2y) + 



B^4+2y(6 + 5y)^ > )e(l-y) 
3 mi mi I 
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- FIG 2 



Figure 2: The electron spectrum for free quark b — > c decay (dashed line), free quark b 
decay (grey line), and B — > X c ev e decay including 1/m 2 corrections (solid line) with Ai = —0.5 
GeV 2 and A 2 = 0.12 GeV 2 . The figure is from '* " 1 



Ai -|- 33A2 r/1 N Ai 



5(l-y)--±6'{l-y) 



3ml 3m 2 



(147) 



Figure || shows the distributions for inclusive semileptonic decays of B mesons. The spec- 
trum close to the endpoint, where the lepton energy becomes maximal, exhibits a sharp spike 
as y — > y m ax- In this region we have we have 



^oce(i-y-p) 



Ai ( p 



(m Q (l-y)) 2 \1 -p 




148) 



which behaves like ^-functions and its derivatives as p — > 0, which can be seen in ( |147| ). This 
behaviour indicates a breakdown of the operator product expansion close to the endpoint, 
since for the spectra the expansion parameter is not 1/mg, but rather l/(m.Q — qv), which 
becomes l/(mg[l — y]) after the integration over the neutrino momentum. In order to obtain 
a description of the endpoint region, one has to perform some resummation of the operator 
product expansion. 

4.5 Resummation in the Endpoint Region 

Very close to the endpoint of the inclusive semileptonic decay spectra only a few resonances 
contribute. In this resonance region one may not expect to have a good description of the 
spectrum using an approach based on parton-hadron duality; here a sum over a few resonances 
will be appropriate. 

In the variable y the size of this resonance region is however of the order of (A/tbq) 2 and 
thus small. In a larger region of the order A/mg, which we shall call the endpoint region, 
many resonances contribute and one may hope to describe the spectrum in this region using 
parton-hadron duality. 



It has been argued in that the £-function-like singularities appearing in ( 147 ) may be 
reinterpreted as the expansion of a non-perturbative function describing the spectrum in the 
endpoint region. Keeping only the singular terms of ( |147|) we write 

~ = 2ft3-2y)S(v), (149) 
T b dy 

where 

00 

S(y) = 9(1 - y) + £ On* (n) (l - y) (150) 

n=0 

is a non-perturbative function given in terms of the moments a n of the spectrum, taken over the 
endpoint region. These moments themselves have an expansion in 1/ttiq such that a n ~ l/mg +1 , 
and we shall consider only the leading term in the expansion of the moments, corresponding to 
the most singular contribution to the endpoint region. 
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Comparing ( 1471) with ( |149|) and ( |150| ) one obtains that 



a 



a i 



dy(S(y)-Q(l-y)) = 0(l/m 2 Q ) 

Ai 



y(S(y) - 9(1 - y)) 



3m2 



(151) 
(152) 



where the integral extends over the endpoint region. 

The non-perturbative function implements a resummation of the most singular terms con- 
tributing to the endpoint and, in the language of deep inelastic scattering, corresponds to the 
leading twist contribution. This resummation has been studied in QCD [|63], |62| and the func- 



tion S(y) may be related to the distribution of the light cone component of the heavy quark 
residual momentum inside the heavy meson. The latter is a fundamental function for inclusive 



heavy-to-light transitions, which has been defined in |62 



i 



2M. 



■{H(v)\h v 6(k+-iD+)h v \H(v)), 



(153) 



where k + = ko + k^ is the positive light cone component of the residual momentum k. The 
relation between the two functions S and / is given by 



S(y) 



l 



dk + f(k^ 



(154) 



from which we infer that the n th moment of the endpoint region is given in terms of the matrix 
element (H(v)\h v (iD + ) n h v \H(v)) . 

The function / is a universal distribution function, which appears in all heavy-to-light 
inclusive decays; another example is the decay B — > X s ^ |64, ^2 , where this function determines 
the photon-energy spectrum in a region of order 1/mg around the K* peak. 

In principle / has to be determined by other methods than the 1/mg expansion, e.g. from 
lattice calculations or from a model, or it has to be determined from experiment by measuring 
the photon spectrum in B — > X s ^ or the lepton spectrum in B —>■ X u lv. In the context of the 



model ACCMM model [51] / has been calculated in [35]. 

Some of the properties of / are known. Its support is — oo < k + < A, it is normalized to 
unity, and its first moment vanishes. Its second moment is given by ai, and its third moment 



has been estimated 



A one-parameter model for / has been suggested in p3, which 



incorporates the known features of / 
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7T 



2 A 



x) exp 



'1 



x 



71 



je(i-x) 



(155) 



where x = k+/A, and the choice A = 570 MeV yields reasonable values for the moments. In 
fig. |3] we show the spectrum for B — > X u £i>£ using the ansatz (|155|) . 

Including the non-perturbative effects yields a reasonably behaved spectrum in the endpoint 
region and the 5-function-like singularities have disappeared. Furthermore, the spectrum now 
extends beyond the parton model endpoint; it is shifted from E™ ax = mq/2 to the physical 
endpoint E™ ax = M#/2, since / is non- vanishing for positive values of k + < A = Mh — uiq. 



30 



FIG 3 



Figure 3: Charged-lepton spectrum in B — > X u £u decays. The solid line is ( |149| ) with the 
ansatz (|155|) , the dashed line shows the prediction of the free-quark decay model. The figure is 
from [B3I. 



5 Conclusions 

The development of the field of heavy quark physics has been indeed remarkable over the last 
few years, experimentally as well as theoretically. From the experimental side, the progress 
in the technology of detectors (e.g. silicon vertex detectors) opened the possibility to study 
b physics even at machines which originally were not designed for this kind of research. In 
this way also the high energy colliders (in particular LEP and TEVATRON) could contribute 
substantially in this area, since they allow to measure states (such as the B s and the b flavoured 
baryons) which lie above the threshold of the T(4s)-5-factories. 

From the theoretical side the heavy quark limit and HQET brought an important success, 
since it provides a model independent and QCD based framework for the description of pro- 
cesses involving heavy quarks. The effective theory approach has originally been formulated 
for exclusive decays but in the past few years a heavy mass expansion has been set up also for 
inclusive transitions. 

As far as exclusive heavy to heavy decays are concerned, the additional symmetries of 
the heavy mass limit restrict the number of non-perturbative functions in a model indepen- 
dent way; furthermore, heavy quark symmetries fix the absolute normalization of some of the 
transition amplitudes at the point of maximum momentum transfer. Phenomenologically this 
has improved our knowledge on the CKM matrix element V c b dramatically; with the value 
\V c f,\ = (39.5 ± 2.0) x 10~ 3 the relative precision of this CKM matrix element is now about 5% 
and thus at a level of the precision with which the Cabbibo angle is known. 

In heavy to light decays heavy quark symmetries do not work as efficiently; in this case only 
the relative normalization of B decays versus the corresponding D decays may be obtained. 
From the experimental side there are first measurements of B — > -niv and B — > piv from the 
CLEO collaboration and an extraction of the CKM matrix element V u b from these processes 
is still to some extent model dependent. The latest value for this CKM matrix element is 
V ub = (3.3 ± 0.2±g;t ± 0.7) x lO" 3 where the last error is due to the model dependence. 

HQET does not yet have much to say about exclusive non-leptonic decays; even for the 
decays B — > D^D^*\ which involves three heavy quarks, heavy quark symmetries are not 
sufficient to yield useful relations between the decay rates pEfl . Of course, with additional 
assumptions such as factorization one can go ahead and relate the non-leptonic decays to 
the semileptonic ones; however, this is a very strong assumption and it is not clear in what 
sense factorization is an approximation. On the other side, the data on the non-leptonic B 
decays support factorization, and first attempts to understand this from QCD and HQET have 
been untertaken f6?fl ; however, the problem of the exclusive non-leptonic decays still needs 
clarification and hopefully the heavy mass expansion will also be useful here. 

The 1/rriQ expansion obtained from the OPE and HQET offers the unique possibility to 
calculate the transition rates for inclusive decays in a QCD based and model independent 
framework. The leading term of this expansion is always the free quark decay, and the first 
non-trivial corrections are in general given by the mean kinetic energy of the heavy quark inside 
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the heavy hadron Ai and the matrix element A2 of the chromomagnetic moment operator. 

The method also allows us to calculate differential distributions, such as the charged lepton 
energy spectrum in inclusive semileptonic decays of heavy hadrons. For this case, the expansion 
parameter is the inverse of the energy release m& — 2E#, where Ee is the lepton energy Close to 
the endpoint, the energy release is small and thus the expansion in its inverse powers becomes 
useless. In this kinematic region one may partially resume the 1/toq expansion, obtaining a 
result closely analogous to the leading twist term in deep inelastic scattering. Particularly in 
the endpoint region a non-perturbative function is needed which corresponds to the parton 
distributions parametrizing the deep inelastic scattering. 

The leading term of the 1/m.Q expansion is the free quark decay, and the result for the 
semileptonic branching fraction in this approximation has been well known for some time. The 
first non-perturbative corrections turn out to be quite small and hence the main corrections 
are the perturbative QCD corrections, where in a recent calculation also the effects of finite 
charm quark mass have been taken into account |68|, |69[| . The radiative corrections lower the 
semileptonic branching fraction somewhat compared to the parton model. 

There has been some discussion on the issue of the semileptonic branching fraction triggered 
by the fact that the data used to be as low as Br{B — > Xtv) ~ 10.5% with a relative error 
of about ten percent. Such a low value for the semileptonic branching fraction in combination 
with tha charm counting in B decays would indicate some theoretical problem; however, the 
recent LEP data yield a value of Br(B — > X£u) = (11.5 ± 0.3)% which is compatible with the 
theoretical expectations. 

The expansion in powers of the inverse quark mass has become the standard tool in heavy 
quark physics and with the forthcomming experiments one may expect a strong improvement 
in our knowledge of the CKM sector of the SM, in particular a test of CP violation as it is 
encoded in the CKM matrix of the SM. 
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Figure Captions 



Figure 1 Diagrams for non-leptonic decays of B mesons. The upper diagram corresponds to the 
leading term in the I/tuq expansion, the middle one to weak anihilation, and the lower 
one to Pauli interference. Diagrams taken from ||56|| . 

Figure 2 The electron spectrum for free quark b — > c decay (dashed line), free quark b — > u decay 
(grey line), and B — > X c ev e decay including 1/m 2 corrections (solid line) with Ai = —0.5 
GeV 2 and A 2 = 0.12 GeV 2 . The figure is from [JTTJ . 



Figure 3 Charged-lepton spectrum in B — > X u iv decays. The solid line is ( |149|) with the ansatz 
( |155| ), the dashed line shows the prediction of the free-quark decay model. The figure is 
from [1631. 
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